Next Article in Journal
Microsyntenic Clusters Reveal Conservation of lncRNAs in Chordates Despite Absence of Sequence Conservation
Previous Article in Journal
Molecular Characterization of Twenty-Five Marine Cyanobacteria Isolated from Coastal Regions of Ireland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Dynamics of the Cell Wall Proteome of Developing Alfalfa Stems

1
Environmental Research and Innovation (ERIN) Department, Luxembourg Institute of Science and Technology (LIST), 4362 Esch/Alzette, Luxembourg
2
Groupe de Recherche en Physiologie végétale (GRPV), Université catholique de Louvain, Earth and Life Institute Agronomy (ELI-A), 1348 Louvain-la-Neuve, Belgium
*
Author to whom correspondence should be addressed.
Considered as equal contributors.
Biology 2019, 8(3), 60; https://doi.org/10.3390/biology8030060
Submission received: 17 June 2019 / Revised: 6 August 2019 / Accepted: 14 August 2019 / Published: 19 August 2019
(This article belongs to the Section Plant Science)

Abstract

:
In this study, the cell-wall-enriched subproteomes at three different heights of alfalfa stems were compared. Since these three heights correspond to different states in stem development, a view on the dynamics of the cell wall proteome during cell maturation is obtained. This study of cell wall protein-enriched fractions forms the basis for a description of the development process of the cell wall and the linking cell wall localized proteins with the evolution of cell wall composition and structure. The sequential extraction of cell wall proteins with CaCl2, EGTA, and LiCl-complemented buffers was combined with a gel-based proteome approach and multivariate analysis. Although the highest similarities were observed between the apical and intermediate stem regions, the proteome patterns are characteristic for each region. Proteins that bind carbohydrates and have proteolytic activity, as well as enzymes involved in glycan remobilization, accumulate in the basal stem region. Beta-amylase and ferritin likewise accumulate more in the basal stem segment. Therefore, remobilization of nutrients appears to be an important process in the oldest stem segment. The intermediate and apical regions are sites of cell wall polymer remodeling, as suggested by the high abundance of proteins involved in the remodeling of the cell wall, such as xyloglucan endoglucosylase, beta-galactosidase, or the BURP-domain containing polygalacturonase non-catalytic subunit. However, the most striking change between the different stem parts is the strong accumulation of a DUF642-conserved domain containing protein in the apical region of the stem, which suggests a particular role of this protein during the early development of stem tissues.

1. Introduction

The plant cell wall is a dynamic structure that modulates a wide variety of physiological events and is a key factor in controlling plant phenotype [1]. This dynamism makes it a barrier against stresses [2], whose properties and structure change during development, depending on the environment [3,4]. The plant cell wall is a composite of polymers (cellulose and heterogeneous hemicelluloses), enzymes, and structural proteins embedded in an aqueous gel of pectins with a variable degree of methyl-esterification. Between the plasma membrane and the primary cell wall of some cells, a secondary cell wall is formed, which, in addition to the carbohydrate-based polymers found in the primary cell wall, may contain variable percentages of non-carbohydrate-based compounds, such as lignin, suberin, or cutin. Furthermore, 5–10% of the total mass of the cell wall is made of proteins [5].
The interest in the study of cell wall proteins (CWPs) has been increasing since the 1990s [5]. Since then, and spurred by the increased use of plant biomass in various industrial fields, the understanding of mechanisms underlying cell wall dynamism and organization has increased [6]. CWPs are mainly secreted to the cell wall after recognition of the signal peptide, which destines them to the endoplasmic reticulum (ER). From the ER, CWPs pass through the Golgi complex and the plasma membrane via secretory vesicles. Although the secretory pathway via the ER-Golgi is the most common, alternative secretory routes have been described [7,8]; a process during which proteins undergo specific maturation events, such as N-linked glycosylations and proline hydroxylation [9,10]. However, recent data indicate that, in addition to the classical signal-peptide-based secretion of proteins to the cell wall, the release of extracellular vesicles to the apoplast may form an important contribution to cell wall remodeling [11]. Restructuration of the cell wall during development is performed by the regulated activity of cell wall-located proteins that target xylans [12], polygalacturonans [13], glycoproteins [14], galacto-oligosaccharides [15], cellulose [16], or lignin synthesis [17]. Moreover, proteins with a domain of unknown function (DUF) were identified in cell wall proteome studies (among others DUF26, DUF231, DUF246, DUF248, DUF288, DUF642, DUF1005, DUF1680), all of which may have a function in cell wall development [16,18].
Proteome studies on plant cell walls have been widely performed on the model species Arabidopsis thaliana. We recently developed a protocol for the extraction of a cell-wall protein-enriched subproteome by combining protocols optimized on A. thaliana and Medicago sativa [19,20]. With this protocol, the contamination with intracellular proteins was limited and a high number of cell wall proteins were extracted [21]. Contrary to A. thaliana, alfalfa is a crop with an economic interest among others, illustrated by the recent creation of a low-lignin variety [22]. We are currently studying this perennial herb as model for cell wall development in response to different abiotic stresses, for instance long-term cadmium exposure [23,24]. Furthermore, we previously studied the differences in gene expression, protein abundance, and cell wall composition between apical, intermediate, and basal alfalfa stem segments [25]. In particular, total proteome analysis revealed a switch in metabolism between the intermediate and the basal segment. In 2012, Verdonk et al. proposed for the first time a comparative cell wall proteome analysis of developmental stages of alfalfa stems [19]. LC-MS/MS analysis indicated that some proteins varied significantly in abundance. These include peroxidase MtPrx 29, GH family 17 (glucan endo-1,3-beta-glucosidase), or polygalacturonase inhibitor [19]. However, the authors state that the quantitative analysis of this study must be interpreted cautiously since it focused on method development and insufficient replicates have been analyzed. Other studies on the cell wall proteome of stems include, among others, studies on A. thaliana [26], sugarcane [27], and flax [28]. Given the current interest in the cell wall proteome, several reviews on the topic have been published, for instance focusing on crops [29]. The abundance of post-translational modifications, such as proline hydroxylation, didehydrophenylalanine, and protein O-glycosylation [10,30,31] and their impact on the function of cell wall localized proteins is a topic of intense research among others for the use of plant systems for protein production.
Here, we describe the experimental work and data from a study aimed at unravelling the significant changes occurring in the cell wall proteome of differently aged tissues, corresponding to different heights along the alfalfa stem. Since the stem of herbaceous plants shows a basipetal lignification gradient accompanying age-dependent progressive secondary growth, the different regions sampled along the stem axis correspond to distinct stages of cell wall maturation. This sampling strategy sheds light on the cell wall dynamics of the growing stem. In summary, our analysis indicates that in the apical region xyloglucan endotransglucosylases, stem 28 kDa glycoprotein, germin-like proteins, DUF642-proteins, and class III peroxidase MtPrx15 are highly abundant. Beta-galactosidases are more abundant in the apical and the intermediate regions. In contrast, endochitinases, aspartic proteinases, peptide-asparagine amidase A (PNGase), class III peroxidases MtPrx29, MtPrx41, and lectins are highly abundant in the mature, basal region. Those proteins that have a significant change in abundance are discussed according to their implication in cell wall metabolism.

2. Materials and Methods

2.1. Plant Material

Soil was collected in a local field (49°33’44′ N, 5°41’49′ E, Musson, Belgium), dried, mixed, sieved at 5 mm, and used to fill 48 containers (7.5 cm × 7.5 cm × 10 cm). Alfalfa (Medicago sativa L. cv Giulia) seeds were inoculated with a peat-based inoculant (HiStick®, Becker Underwood), according to the manufacturer’s instructions, and sown 3 seeds per pot. Plants were grown in incubators (Fitotron SGC 120, Weiss Technik UK, Leicestershire, UK) at 22 °C/17 °C, 13 h/11 h d/n, 60% humidity. Following germination, 2 of the 3 plants per container were selected and watered with deionized water. A first cut was performed 2 months after sowing, and plants were allowed to regrow. Soils were then weekly fertilized with 50 mL macronutrient solution (Hoagland 1X macronutrients only) [32]. All growing shoots reaching 25 cm and developing at least 9 internodes were sampled 6 weeks after the first cut. Stems were separated from leaves and stipules with a razor blade and divided into 3 parts of equal length, namely the apical, intermediate, and basal regions. In a previous study [25] the stem was divided into five parts and part 1/5, 3/5, and 5/5 were analyzed to have more discriminative samples. The requirement of larger samples for cell wall protein enrichment necessitated the application of a different sampling strategy. For each stem segment, 4 individual biological replicates (4 g to 5 g for fresh weight) were prepared and stored at −80 °C prior to protein extraction.

2.2. Cell Wall Protein Enrichment

2.2.1. Extraction of Cell Wall Proteins

Cell wall enrichment was attained using a sucrose gradient procedure as described by Feiz et al. [20]. Protein extraction was done with a sequential extraction method using buffers supplemented with CaCl2, ethylene glycol-bis(β-aminoethyl ether)-N,N,N′,N′-tetraacetic acid (EGTA), and LiCl [21].
The frozen plant material was ground in liquid nitrogen and transferred to a 50 mL Falcon® tube with 20 mL 0.4 M sucrose, 5 mM Na+ acetate, pH 4.6 at 4 °C. After vigorous shaking (2 min at 25 Hz), the samples were placed on a rocking platform (overnight, 4 °C) and subsequently centrifuged (1000 g, 15 min, 4 °C). After centrifugation, the supernatants were discarded and the pellets were resuspended in 10 mL of the same buffer, but with 0.6 M instead of 0.4 M sucrose. After 30 min on a rocking platform at 4 °C, the samples were centrifuged (1000 g, 15 min, 4 °C) and the supernatants were again discarded. This washing step was repeated respectively with the same buffer with 1 M sucrose and, subsequently, twice with the buffer without sucrose. The pellets, corresponding to the cell wall-enriched samples, were transferred to 30 mL tubes.
To these cell wall-enriched samples, 7.5 mL 5 mM Na+ acetate, 200 mM CaCl2, and pH 4.6, 4 °C was added, and the samples shook on a rocking platform for 30 min at 4 °C. After centrifugation (10,000 g, 15 min, 4 °C), the supernatants were saved, the procedure repeated, and supernatants pooled. These pooled supernatants formed the CaCl2 fraction.
A total of 10 mL 5 mM Na+ acetate and 50 mM EGTA (pH 4.6) was added to the pellet, which shook vigorously for 1 h at 37 °C. After centrifugation (10,000 g, 15 min, 4 °C), the supernatants were saved and the incubation with EGTA was repeated twice. The EGTA fraction is composed of the pooled supernatants of these three repetitions.
A total of 15 mL 5 mM Na+ acetate and 3 M LiCl, (pH 4.6) was added to the pellet after EGTA extraction, The pellet was resuspended and placed overnight on a rocking platform (4 °C). The supernatants after centrifugation (10,000 g, 15 min, 4 °C) were saved, forming the LiCl fraction.

2.2.2. Concentration and Desalting of the Extracts.

Amicon Ultra-15 10K Centrifugal Filter Devices (Millipore) were used to reduce the sample volume to approximately 200 µL (4700 g, 4 °C). This final volume was washed and desalted using the Bio-Rad ReadyPrep 2-D Cleanup Kit. Finally, the dried protein samples were solubilized in 100 µL labelling buffer (7 M urea, 2 M thiourea, 2% (w/v) CHAPS, 30 mM Tris), and the Bradford protein assay with bovine serum albumin as standard was used to determine the protein concentrations.

2.3. Cell Wall Protein Analysis by 2D-DIGE

2.3.1. Protein Migration

Distinct gel images were obtained for the three factions. Therefore, these were analyzed separately, using a separate internal standard and separate analysis. Labelling, migration, and spot detection were performed as previously described [25] with the following slight modifications: (1) 50 μg of proteins for the CaCl2-fractions, 40 μg of proteins for the EGTA, and 40 μg of proteins for the LiCl fractions were labelled with respectively 400 pmol and 320 pmol of the cyanine NHS-dyes LUMIPROBE LLC (Hannover, Germany); (2) following labelling, resolubilization in the lysis buffer and an addition of Destreak Reagent (GE-Healthcare, Chiago, IL, USA) and Ampholyte (Bio-Rad, Hercules, CA, United States) (pH 3–10) samples were loaded onto ImmobilineTM DryStrip (3–10 NL, 24 cm) for overnight rehydration; (3) a 5 step program was used for isoelectric focusing (IEF): (1) 100 V for 3h, (2) a 4 h linear gradient from 100 V to 1000 V, (3) 1000 V for 6 h, (4) linear gradient from 1000 V to 10,000 V for 6 h, and (5) 10,000 V till 95,000 Vh were reached. During IEF, the current was limited to 75 μA per strip. The use of different protein quantities for the different fractions has no influence on normalization and quantifications since the three fractions were analyzed independently.

2.3.2. Spot Selection

Following spot detection and matching using Decyder (v7.0, GE-Healthcare), the normal volume of each spot was extracted, and statistical analysis was performed using R (v3.1.2) with the package FactoMineR v1.0. Prior to statistical analysis, the normal spot volumes were standardized by dividing the volume of each spot by the corresponding spot volume of the internal standard. Spot average ratio r (Apex vs Intermediate, Apex vs Base, and Intermediate vs Base) was calculated. In case this r ratio was inferior to 1, the additive inverse of the inverted r (−1/r) is reported. Data (standardized spot volumes) were further log2 transformed. Only spots with no more than 1 missing value per condition were kept for analysis. Comparisons (Apex vs Intermediate, Apex vs Base, and Intermediate vs Base) were performed by spot, using a two-tailed Student t-test for samples of equal variance. Only spots with a r-value <−1.5 or >1.5 and p-values inferior to 0.05 were considered to vary significantly. Spots of interest were visually checked using the Decyder (v7.0, GE-Healthcare) spot view and selected for picking.

2.3.3. Spot Picking, Digestion, MS/MS Analysis

Selected spots were picked with an Ettan Spot Picker (GE- Healthcare). Digestion was carried out using a Freedom EVO II workstation (Tecan, Männedorf, Switzerland). Briefly, gel plugs were washed with 50 mM ammonium bicarbonate in 50% (v/v) methanol/MQ water and dehydrated with 75% acetonitrile (ACN). After dehydration, proteins were digested with trypsin Gold (Promega, Fitchburg, WI, USA) and 8 µL of a solution containing 5 ng/µL trypsin in 20 mM ammonium bicarbonate (overnight, 37 °C). After digestion, peptides were extracted from the gel plugs with 50% (v/v) ACN/0.1% (v/v) trifluoroacetic acid and dried. Peptides were solubilized in 0.7 µL of 50% (v/v) ACN/ 0.1% (v/v) TFA and spotted on MALDI target plates. To this, 0.7 μL matrix (7 mg/mL α-cyano-4-hydroxycinnamic acid in 50% (v/v) ACN/0.1% (v/v) TFA) was added.
A peptide mass spectrum was acquired using the 5800 MALDI-TOF/TOF MS (Sciex, Darmstadt, Germany), the 10 most abundant peaks, excluding known contaminants, automatically fragmented. MS analyses were carried out as described by Printz et al. [33]. MS and MS/MS spectra were submitted for database-dependent identification against the NCBInr database limited to the taxonomy Viridiplantae (1,717,798 sequences) (http://www.ncbi.nlm.nih.gov) on an in-house MASCOT server. A second search was done against an EST Fabaceae database (19,932,450 sequences). The parameters used for these searches were mass tolerance MS 100 ppm, mass tolerance MS/MS 0.5 Da, maximum 2 missed cleavages, fixed modification carbamidomethyl-cysteine, variable modifications oxidation of methionine, oxidation of tryptophan kynurenine, and double oxidation of tryptophan. Proteins were considered identified when at least two peptides passed the MASCOT-calculated 0.05 threshold scores (Supplementary Table S1A–G). All identifications were manually validated. When high quality MS/MS spectra were not identified in database searches, the sequence of the peptides was determined manually. In order to obtain an objective confidence score for these peptides, the spectra were resubmitted with adjusted search parameters. After manual validation, these were added in the Supplementary Table S1A–C with the peptide score in red. When peptides in the same spot matched different database entries, they were aligned, and it was verified that these belong to the same protein using the NCBI BLAST and the Alfalfa Gene Index and Expression Atlas Database (AGED) databases (http://plantgrn.noble.org/AGED/). When a protein with a trivial name was identified, a BLAST alignment against the viridiplantae database was done and the name of the protein with the highest significant homology was added in the result table. For the prediction of the localization of the proteins, the recently developed algorithm DeepLoc was used [34].

2.3.4. Post-Identification Statistical Procedure

Spots containing 1 significant protein (116) were analyzed using the freeware R (v3.1.2) with the additional package FactoMineR v1.0. Principal component analysis (PCA) was performed on this set of proteins (Supplementary Table S2A). PCA was carried out with the R-function PCA by using the normalized spot volumes of the 116 spots in each stem region. Proteins were clustered in 10 groups with the hierarchical clustering on principle components (HCPC). The R function and PCA axis were characterized using the dimdesc function of the FactoMineR v1.0 package [35]. This procedure allowed the extraction of a set of 103 proteins that highly contribute (p < 0.05) to Dim1 and/or Dim2, among which 61 (61.6%) and 4 (66.7%) proteins were predicted to be secreted (Supplementary Table S2B,C). Peroxidases were identified using the web-tool PeroxiBase, available at http://peroxibase.toulouse.inra.fr/ [36]. Multiple sequence alignment was performed using Clustal Omega, available at: http://www.ebi.ac.uk/Tools/msa/clustalo/.

3. Results

The cell wall-enriched proteome of the stem of alfalfa (1st year of growth, 2nd cut from plants with 3 to 4 shoots) was extracted from three different stem regions: The apical, intermediate, and basal stem. CWPs were extracted according to a protocol optimized on alfalfa stems and derived from previous work on plant cell wall proteomics [19,20,21]. This procedure includes a first step of cell wall enrichment followed by a 3-step sequential extraction procedure with buffers containing CaCl2, EGTA, and LiCl. Each of these fractions was analyzed separately, and normalization and the following quantification was done for the CaCl2, EGTA, and LiCl fractions independently. Four biological replicates were generated for extraction. Proteins were labelled using fluorescent dyes and separated on 2D-gels (2D-DIGE) prior to quantitative analysis and protein identification (Supplementary Figure S1), as previously done by Watson et al [37]. Since no current proteomics techniques allow the analysis of all proteins present in a sample, the here-reported results correspond to that part of the proteome that can be analyzed with the used approach and is limited by the extraction protocol and the use of 2D-DIGE for visualization and quantification. The intracellular proteins identified are reproducibly found and their abundances changed significantly. They are part of the subproteome that is extracted and thus are discussed as such. The big groups of proteins identified in this study are similar to those found in other cell wall proteome studies [27,38,39].
Only spots matched on 75% of the gels and displaying a volume ratio r superior to 1.5 or inferior to −1.5 (Student t-test p-value <0.05) were considered as significant. According to these criteria, 91 significant spots were selected and analyzed in the CaCl2 fraction, 139 in the EGTA fraction, and 27 in the LiCl fraction. The MS analysis of these 257 spots resulted in 186 significant identifications, among which 122 (65.6%) were predicted to carry a signal peptide (SignalP) [40]. Since the presence of more than one protein in a spot excludes the quantification of the individual protein, only the 116 spots, wherein exactly one protein was identified, can be used for biological interpretation. Given that a different location on the gel (change in pI and/or molecular weight) of two spots, nominally containing the same protein, indicates that these are either isoforms encoded by different, homologous genes or different forms of the same gene product due to post-translational processes or artefactual changes, which are treated separately. The use of DeepLoc on these 116 spots resulted in the prediction of secreted, soluble protein for 65 (56%) spots. Although Deeploc performs well in predicting the subcellular localization of proteins, only the actual observation of a protein outside of the plasma membrane can confirm this. The change of one amino acid or, for instance, saturation of the vacuolar sorting system [41] are among the factors that may influence the actual localization of a protein.
These cell wall-predicted proteins were classified based on previously listed functional classes [42]. This classification was performed according to the list of domain hits proposed by the NCBI blast web-tool and from available literature (Supplementary Table S2D). The most represented class of proteins that undergo significant variations in abundance along the stem of alfalfa is “proteins acting on carbohydrates”. The other proteins are essentially classified as “oxido-reductases”, “proteases”, “proteins with interaction domains”, “proteins with nutrient reserve function”, and “miscellaneous proteins”. The highest fold-change of the study (7.56, Apex vs Base) was found for an unknown protein (NCBInr gi:537313) 98% identical in sequence with the Medicago truncatala plant F18G18-200 protein (CM001217.2) with a DUF642 conserved domain that was identified in four spots (two in the CaCl2 fraction (1586 and 1609) and two in the EGTA fraction (956 and 965)). In the following discussion and tables, the name F18G18-200 is used for this protein. The protein identified in spot CaCl-2827 is 90% identical in sequence to acidic endochitinase and discussed as such.
PCA followed by post-PCA clustering in 10 groups was performed on the set of 116 spots, in which only one protein was identified (Figure 1 and Table 1). The first component (Dim 1) of the PCA explains more than half (55.98%) of the variability of the cloud of points and allows us to differentiate the three stem regions, indicating that each stem segment has a characteristic cell wall proteome. Similar discrimination was obtained in a non-targeted proteome study of the same three stem regions [25]. The highest similarities were observed between the apical and intermediate stem regions. The second dimension of the PCA accounts for 10.24% of the variability and differentiates the intermediate from the basal and apical stem regions. However, this dimension is represented by only six spots, which significantly contribute to this axis (p < 0.05), indicating that few proteins had a behavior specific to this segment. From the 116 proteins that changed significantly, a subset of 103 proteins was identified as “high contributors” to these two PCA-axes. These were extracted using the dimdesc function of ‘FactoMineR’ at the threshold of 0.05 (Supplementary Table S1B,C). From this selection, two distinct trends emerged. The proteins either accumulated more in the apical (49 proteins) or basal (50 proteins) part of the stem. Only four proteins, all from the cysteine protease family (CaCl2 2269, EGTA 1450, EGTA 1452, and EGTA 1451) accumulated mainly in the intermediate stem region.
Proteins showing a higher accumulation in the most lignified basal region were grouped mainly in clusters 1 to 4 (Table 1) and were identified, among others, as members of the lectin family: PNGase, chitinases/endochitinases, and glucan endo-1,3-β-glucosidases (Table 1 and Supplementary Table S2D). Proteins not predicted to be cell wall localized also grouped in these clusters, especially ferritins, beta-amylases, and glutelins type-A (Table 1 and Supplementary Table S2D). In contrast, proteins belonging to groups 6 to 10 were more abundant in the young, apical segment. The most remarkable changes in comparison with the base of the stem were the high abundance of beta-galactosidases, germin-like proteins, class III peroxidases, plant/F18G18-200 proteins, stem 28 kDa glycoproteins, polygalacturonase non-catalytic proteins, and xyloglucan endotransglucosylases. Confirming findings that chloroplast-located proteins are more abundant in the green, apical region, it was found that non-cell wall proteins from these clusters are chloroplast localized. We further noted the higher abundance of one dirigent-like protein (Pfam 03018), putatively involved in stereoselective lignan biosynthesis (Table 1).
The abundance of these 103 high contributor proteins in the intermediate stem region was generally intermediate between that in the two extreme stem parts. However, proteins from cluster 7, which contained 18 isoforms of beta-galactosidase, accumulated to a similar extent in the apical and intermediate region. Proteins from cluster 1, containing mainly lectin-like proteins, accumulated in the basal segment and had a lower abundance in the younger stem parts. In contrast, the two DUF642 proteins (Plant/F18G18 in Table 1) from cluster 10 strongly accumulated in the apical part. Figure 2 summarizes the most important changes occurring along the alfalfa stem identified in this study.
During MS analysis, several post-translational modifications were observed and are added in the Supplementary Table S1A–C, most of which concern the confirmation of predicted signal peptides. The proteolytic cleavage of N-terminal presequences or inhibitor domains (for instance the spots 914, 969, and 1233 in the EGTA fraction) was confirmed by a literature search or blasting the found sequence and identification of the active domain. Finally, α-β dehydration of phenylalanines in the beta-subunit of polygalacturonase (LiCl spot 300), as recently described [30], was confirmed.
A C-terminally truncated calreticulin was observed (CaCl2 spot 1274). Although regarded as an ER-resident protein, calreticulin was described in other subcellular compartments [43]. Animal calreticulin is linked with the recognition of tumor and apoptotic cells, and wound healing in plants it is linked to growth and stress responses [44]. Although a mass loss of approximately 6 kDa was observed in vitro [45], to our knowledge, cleavage in the C-domain of calreticulin was not yet observed in vivo. The found cleavage site DPVD is conserved in plant homologues and corresponds to a caspase-like activity, a regulatory proteolytic activity that is known in plants [46]. Although our data explains the lack of immunological detection of the C-terminal HDEL-sequence for plasma membrane-linked calreticulin [47], a further study should clarify whether this observation stems from processing or degradation.

4. Discussion

Uniform stems of alfalfa were dissected in equal thirds and cell wall protein-enriched fractions extracted using a sequential extraction protocol. The different sub-proteomes were independently analyzed using a gel-based approach. This quantification strategy was selected because the fractions contain relatively few different proteins, thus limiting the masking effect when 2D gels are used on more complex samples. Furthermore, a separation of proteins, and not peptides as in LC-MS-based approaches, allows visualizing the heterogeneity of the different forms, in which a gene product can be present in a sample. This has the drawback of having more than one spot and thus a quantitative value for proteins that are heterogeneous, hindering straightforward interpretation. However, this heterogeneity is biological. A way to overcome the inherent limits of gel-based proteomics while keeping the capacity to visualize all different forms in which a protein is present is the use of top-down proteomics, an approach not yet ready for application on the system studied here [48].

4.1. Carbohydrate-Interacting Proteins Accumulate in the Cell Wall of the Apical Region

The cell wall of the apical region is characterized by a significantly higher presence of a stem 28 kDa glycoprotein (EGTA-1391), a polygalacturonase non-catalytic protein (LiCl-300), and xyloglucan endoglucosyltransferase (Glycoside Hydrolase 16 (GH16)—CaCl2-1945 and CaCl2-1940), which cluster in groups 8 and 9 (Table 1). In the alfalfa stem, we previously observed that the abundance of this 28 kDa glycoprotein varied significantly with changes in copper availability in the substrate [49], and Verdonk et al. also observed a higher accumulation of this protein in the apical region [19]. In particular, it accumulated more in the apical region when copper availability was optimal. Here, the protein is found in the highest abundance in the youngest stem region (containing elongating internodes), which correlates with a high expression level of the contig 59881 (plantgrn.noble.org/AGED). This corroborates the attribution of a specific role to this protein during the early growth of the stem [50].
The higher abundance of the polygalacturonase (PG) non-catalytic protein (β-subunit) (LiCl-300) is consistent with the higher expression of contigs coding for this protein. The data on the contigs 111493, 12933, 27771, 48302, 51002, 22692, 53836, 12932, 12931, and 65583—all coding for isoforms of this protein—as can be found at http://plantgrn.noble.org/AGED/, indicate a higher expression in younger tissues. The expression data thus reinforces our data, indicating that these proteins are involved in early stem development. The function of these proteins remains elusive. Tomato plants expressing a PG β-subunit antisense transgene show an increased degradation of the pectin fraction of the cell wall [51]. It was thus hypothesized that the presence of the β-subunit restricts pectin solubilization and limits pectin depolymerization during fruit ripening. In particular, the β-subunit may protect or limit the access of the catalytic PG protein to pectin sites within the cell wall. The hypothesis that PG β-subunit modulates in vivo PG activity was supported by findings on wound-inducible PG β-subunit antisense tomato plants, which resulted in elevated PG activity [52]. In the young tissues, such as the elongating stem of alfalfa, the high abundance of PG non-catalytic proteins may thus limit pectin degradation.
The two xyloglucan endotransglycosylase/hydrolase (XTH) (CaCl2-1945 and CaCl2-1940) share respectively 80% and 79% identity with A. thaliana protein AT2G06850 and 86% with VrXTH1 from Vigna radiata. In the A. thaliana uro mutant, with delayed secondary cell wall formation, the expression of XTH-coding genes (including AT2G06850) is significantly reduced [53]. In V. radiata, the VrXTH1 gene is inducible by brassinosteroids, auxin, and Ca2+ and is highly expressed in the hypocotyl and stem, more specifically in the epidermis and cortical layers, where cell elongation occurs [12,54]. XTH4 catalyzes the hydrolysis of ß(1–4) bonds in the xyloglucan backbone and transfers the xyloglucanyl segment to the non-reducing glucose residue of an acceptor [55], being a xyloglucan or a xyloglucan oligosaccharide. This protein may thus ease cell growth by modifying the cell wall structure through cleavage and rejoining of cross-links between xyloglucan and microfibrils. The high abundance of these proteins in the apical stem region is in agreement with observations done in V. radiata hypocotyls [12].
Young stem tissues are also characterized by the high abundance of DUF642-domain proteins, with an unknown function (NCBInr gi:537313, homologous to Medtr1g011800.1) [56]. In 2016, it was mentioned that proteins with this domain were found in all cell wall proteome studies [27]. These proteins were identified in both the CaCl2 (spots 1586 and 1609) and EGTA (spots 956 and 965) fraction and underwent the most important change in accumulation in this study, indicating a role of DUF642-domain proteins in the development of younger stem segments. In A. thaliana, BIIDXI, a DUF642 protein, is linked with hypocotyl growth via the efflux of auxin [57]. Whether other DUF642 proteins have similar functions remains to be determined. In a study on the development of etiolated A. thaliana hypocotyls, few changes were observed in the accumulation of five proteins with a DUF642-domain [38]. For at least some isoforms of DUF642-proteins, in vitro interaction with cellulose and, to a lower extent, hemicellulose, has been observed [16]. Phylogenetic analyses have shown that this domain is specific to the angiosperm and gymnosperm lineage [16]. The differential expression of DUF642 genes in different plant organs and at different time points indicates the existence of a tight spatiotemporal regulation of these genes and suggests that different members may exert a distinct role in plant development [58,59].
Doing an in silico analysis, we identified 7 M. truncatula genes coding for homologous DUF642 proteins (Medtr1g011800.1, Medtr2g103170.1, Medtr0045s0110.1, Medtr4g039680.1, Medtr4g039720.1, Medtr4g039740.1, Medtr2g019600.1). All seven genes have one to six GTAC-motifs in their 1000 bps 5’-UTR promoter region, a motif which may be recognized by the SQUAMOSA promoter-binding protein-like 7 (SPL7) transcription factor, the central regulator of Cu-homeostasis in plants [60]. In addition to that, a promotor analysis of these seven genes revealed two conserved motifs (Figure 3). Using the JASPAR CORE 2016 plants database resulted in the identification of EDT1/HDG11 and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) as candidates for binding to these promoter-motifs. EDT1/HDG11 is a HD-ZIP IV transcription factor upregulating cell wall-loosening genes [61]. Therefore, DUF642 may be regulated via the same HDG11-dependent transcriptional circuit and may partake in elongation in young stem regions. (SOC1) is a MADS-box transcription factor involved in flowering time control. When induced in the shoot apex, it regulates (via gibberellin interaction and photoperiod perception) the transition from vegetative to reproductive growth. In this respect, it should be noted that SPL3 is directly regulated by SOC1 [62]. Furthermore, some members of the DUF642 protein family directly interact with pectin methylesterase (PME) isoforms, generating the hypothesis that they modulate the pectin methylesterified state in vivo and are thus involved in the fine-tuning of the biomechanical properties of the cell wall [18,58]. The involvement of DUF642 proteins in cell elongation in young tissues was recently shown [63], and the identification of these conserved promoter motifs as binding sites for transcription factors involved in cell wall loosening may contribute to the further disentanglement of their physiological function.

4.2. β-galactosidases are Abundant in the Apical and Intermediate Stem Part

Two β-galactosidase proteins (NCBInr gi:657400954 (one spot) and NCBInr gi:357437609 (18 spots), Glycoside Hydrolase 35), sharing respectively 70% and 67% identity with the A. thaliana homologue AT4G26140.1 (Gal-12, β-s subfamily a1), accumulate more in the intermediate and apical region compared to the basal segment. A similar increased accumulation of beta-galactosidase was previously observed in the younger tissues of alfalfa stems compared to older internodes [19]. Moreover, the higher expression of genes coding for β-galactosidases in young tissues was previously observed. For instance, Iglesias et al. (2006) report a higher expression of AT4G26140 in young leaves and apical stems compared to mature leaves and the basal stem region [64]. Likewise, a transcriptome study on alfalfa found a more than five-fold higher expression of beta-galactosidase genes in younger compared to older tissues [65]. As all members of this family, Gal-12 hydrolyses β-(1,4)- and β-(1,3)-linked galacto-oligosaccharides, such as those found in the pectins galactobiose and galactotriose. However, since this protein is unable to cleave galactose from xyloglucan oligosaccharides, larch wood arabinogalactan, or oat xylan, it is unlikely to be implicated in hemicellulose remodeling [15]. In our study, high accumulation of β-galactosidases is observed in the apical and intermediate stem region, confirming its implication in the youngest regions. Since mature root sections labelled with Gal-12 antibodies showed specific labelling in thickened walls of xylem cells, its homologue, which accumulates in young alfalfa stem regions, may act in the remodeling of the xylem cell wall.

4.3. Class III Peroxidases MtPrx29, MtPrx41, and MtPrx15

Whereas most CWPs from the same family follow similar variations in abundance, different members of the class III peroxidase family (Class III Pfam PF00141) behave differentially. Similar disperse changes in abundance during cell wall maturation were previously observed [38]. Proteins identified with the two NCBInr identifiers gi:357491415 (MtPrx29/MsPrx05—cluster 5, CaCl2-1769) and gi:657377089 (MtPrx41—cluster 4, EGTA-1168) accumulate more in the basal region, an observation previously made when developing alfalfa [19]. The corresponding contigs (contig 65262, 88968, and 16417) have a higher expression in post-elongating internodes. In contrast, the three spots wherein the protein was identified with the NCBInr identifier gi:537317 (MtPrx15/ MsPrx15—cluster 8, LiCl-340, LiCl-1156, LiCl-1163) are of higher abundance in the apical region, and since they are identified in the LiCl extract, they interact more strongly with the cell wall. Class III peroxidases belong to a multigenic protein family, involved in processes such as limiting the spread of infections to the strengthening of the cell wall [66]. During plant growth, cell expansion relies on a balance between cell wall stiffening and cell wall loosening, antagonistic processes wherein members of the class III peroxidase family are involved [67]. Although CW peroxidases belong to large multigenic families (139 members in Medicago truncatula [37]), only a few members have been functionally characterized in planta [67], a task made more difficult by their functional redundancy. Targeted gene expression analysis under different growth conditions provides an alternative strategy to study their role in planta. In this respect, in a recent study we identified some alfalfa class III peroxidases whose expression is affected by abiotic stresses [68].
The three MsPrx15 isoforms we identified share the highest similarity with the A. thaliana class III peroxidase AtPrx53. However, given the limited sequence identity of 58%, no specific function can be attributed. Unfortunately, data concerning MsPrx15, or its closest homologue MtPrx15, are scarce. In physiological conditions, a contig showing high sequence similarity to MsPrx15 (contig_12706, e-147) is upregulated in elongating internodes. Rahman et al. (2015) showed that this apoplast-located peroxidase accumulates in roots of the salt-tolerant alfalfa ecotype NM-801 following salt treatment [69]. In Lotus japonicus, inoculation with Mesorhizobium loti induced the down-regulation of the gene coding for a MtPrx15-homologue in the roots of the plant. In M. truncatula, leaf wounding induces a reduction in the abundance of MtPrx15, which can be recovered by application (3 min after wounding) with the redox inhibitor diphenyleneiodonium, suggesting that MtPrx15 abundance varies with the presence of reactive oxygen species [70]. This peroxidase is furthermore a target of cysteine sulfenylation, indicating that this protein may be redox-regulated [71].
Recently, sequence variations and glycosylation events occurring on three different isoforms of this peroxidase were identified [25], without identifying changes that explain the pI-shifts observed on 2D-gels. This peroxidase is the major protein present in the LiCl fraction of alfalfa stems (representing half of the total volume of all spots on the gels of the LiCl fraction). In the current study, the higher accumulation of MtPrx15-isoforms in the upper stem regions, a site of broad metabolic activity directed toward cell wall development, indicates a role in the development of the cell wall, probably by promoting the cross-linking of lignin precursors.
A dirigent-like protein (CaCl2-2570) also attains its highest accumulation in the youngest stem region. The expression of the closest alfalfa homologues of the identified protein (AGED database contig_5031 and contig_29341) are likewise higher in young, elongating stems compared to post-elongation stems. A concerted action of MtPrx15 and dirigent-like protein in lignin biosynthesis or the synthesis of lignans can be hypothesized. Overexpression of a soybean homologue resulted in an increased accumulation of lignans in Pisum sojae and a concomitant increase pathogen resistance [72]. The shunt of monolignols to lignans via the intervention of dirigent-like proteins may therefore increase pathogen resistance, as is found in other species [73,74]. In contrast, the two class III peroxidases MtPrx29 and MtPrx41, which accumulate more in the oldest stem regions, may exert a different role in cell wall metabolism, but this latter has not been yet elucidated.

4.4. Enzymes Involved in Protein Degradation/Maturation Accumulate in the Oldest Stem Region

With increasing tissue maturity, a higher accumulation of a PNGase (CaCl2-1013) (peptide-asparagine amidase A, glycopeptidase, N-oligosaccharide glycopeptidase) is observed. In a cell wall proteome study on Brachypodium distachyon grains, a similar higher abundance of PNGase was observed in more mature cell walls [75]. This protein hydrolyses N-glycans with α-1,3-fucosylated asparagine-bound GlcNAc, as commonly found in secreted cell wall proteins [14,76]. At least two forms of PNGases can be discerned. Cytosolic PNGase (cPNGase) hydrolyses amide linkages in glycosylated Asn-residues of misfolded proteins prior to proteosomal degradation and have an optimum activity in the neutral pH region. The plant-specific acidic PNGase (aPNGase) identified in this study has an acidic optimum pH and hydrolyses glycopeptides carrying high-mannose type or complex type N-glycans in the vacuole and cell wall [77]. The aPNGase (NCBInr gi:922395483) was found in the CaCl2 fraction and decreased in abundance from the lignified basal region to the younger stem segments. Alignment of this aPNGase with the one from almond and tomato (referred to as PNGase-A and PNGase-Le) revealed that the proteins share three conserved regions (Supplementary Figure S2). These may be important in determining the activity of this protein. By catalyzing the hydrolysis of N-glycans from glycoproteins, its accumulation in the oldest stem region suggests that free N-glycans are preferentially generated in the apoplasm of the oldest tissues [77]. A role in the recovery of sugars from proteins in aging tissues can thus be considered.
This hypothesis would corroborate with the higher accumulation of different endo- and exoproteases in the basal region. These include aspartic proteinase Asp1 (EGTA-1119), serine carboxypeptidase S10 (EGTA-1111), eukaryotic aspartyl protease family protein (EGTA-914), and cysteine proteinase family members (EGTA-1228, EGTA-1233, and CaCl2-2828). However, other cysteine-proteinase family members (including vignain-like proteins) behave differently (CaCl2-2269, EGTA-1296, EGTA-1450, EGTA-1451, EGTA-1452), indicating the differential regulation of members of this protein family according to the developmental maturity. A functional characterization of the different proteases identified here is currently not available. By the careful examination of the mass spectra, we were able to show that the N-terminal inhibitor domain of the eukaryotic aspartyl protease gi:657371151 (EGTA-914) and the cysteine proteinase gi:357437719 (EGTA-1228, EGTA-1233) is cleaved off (Supplementary Table S1B). The protein identified in the spots EGTA-1228 and EGTA-1233 is homologous to cysteine proteinase RD21A, a protein previously identified in the apoplastic fluid of A. thaliana [78], in seed mucilage [79], and in the culture medium, wherein A. thaliana seedlings have grown [80]. This defense-related protein is predicted to undergo a three-step maturation; cleavage of the signal peptide, cleavage of N-terminal inhibitory domain (observed as the peptide with m/z 1376 Supplementary Table S1B), and finally cleavage of the C-terminal granulin domain [81]. Although we did not identify the peptide corresponding to the latter cleavage, the position on gel (Supplementary Figure S1) and the lack of peptides covering the granulin-domain while the protease domain is well covered suggest that the protein is in its fully active form. Interestingly, RD21A homologous from other species have been found in the apoplast [82], and recently the retargeting of this protein under the influence of the cyst nematode effector 4E02 was described [83].
The abundance of chloroplastic ferritin 2 and 3 was significantly higher in the basal and intermediate regions. Ferritins are ubiquitous proteins that form complexes which can gather up to 4500 atoms of iron [84] and are of major importance for iron homeostasis. In 2009, Ravet et al. showed that ferritins mainly function as defense mechanisms against free iron-induced oxidative stress and not as storage for iron needed for plant development [85]. In Lupinus luteus, the accumulation of ferritin mRNA and the corresponding protein was reported in senescing nodules, supporting the idea that this protein participates in iron translocation from senescing tissues to other plant regions [86]. The presence of four (NCBInr gi:357468557) and two (NCBInr gi:357492793) GTAC motifs in the 1000 bps 5’-UTR sequence upstream of the genes coding for these two2 proteins in M. truncatula indicates that ferritin expression may be regulated by the SPL7 transcription factor.
Remobilization of metabolites from the mature stem region was further suggested by the higher accumulation of nine beta-amylase isoforms (NCBInr gi:3913031 (homologous with Medtr4g045667.1)), involved in starch remobilization [87] and the basal regions, in comparison with both the intermediate and the apical segments. Interestingly, the screening of the 1000 bps 5’-UTR sequence upstream of the gene coding for this protein (Medtr4g045667.1) revealed the presence of four GTAC motifs, making it also potentially dependent on SPL7 regulation. These observations suggest that copper availability affects the accumulation of proteins typically present at certain stages of stem maturity. Consequently, stem development might be partially dependent on copper availability.
We finally noted the higher abundance of two LecRK-domain containing lectins (respectively NCBInr gi:2951684 (CaCl2-2127, CaCl2-2128, EGTA-1411, LiCl-765) and gi:400180 (CaCl2-2193, EGTA-1390, EGTA-1420) in the cell wall fraction for the basal stem. Lectin-domain proteins are known to be involved in the defense against pathogens and herbivores [88], and ATP-binding LecRK homologues are regarded as having antimicrobial activities [89].

5. Conclusions

This study reveals the impact of tissue maturity on the cell wall protein-enriched subproteome extracted by a sequential extraction protocol. Stem sections at three different heights along the stem, and thus different stages of cell wall development, were studied. The strongest similarities are observed between the apical and intermediate stem sections, both of which contain a high abundance of proteins that target and modulate the structure of cell wall polymers. DUF642-containing proteins, beta-galactosidase, the polygalacturonase β-subunit, stem 28 kDa glycoproteins, and endoxyloglucan transferase appear to have an important function in the youngest regions, confirming previous data [19]. These proteins are involved in modifying carbohydrate-based polymeric backbones, therefore, their higher abundance relates to the transition from young elongating regions to mature, thickening tissues.
In the most basal, oldest stem region, proteins such as PNGase and beta-amylase accumulate, suggesting that in more mature stem tissues the cell wall proteome is more oriented towards the recycling of sugars from glycoproteins and storage molecules. By studying the cytosolic proteome of the stem [25], we previously found that defence-related proteins accumulate in the basal region. Here, we report in addition the accumulation of a defensin precursor protein, various chitinases/endochitinases, and members of the class III peroxidase-family in the most mature stem segment. These observations confirm that protein-based defence is important in mature alfalfa stem cell walls. Opposed to this, the higher accumulation of dirigent-like protein in the younger tissues suggests that defence-activities are performed by the production of lignans.

Supplementary Materials

The following are available online at https://www.mdpi.com/2079-7737/8/3/60/s1. Table S1A–G. Protein identifications and associated mass spectrometric data. Table S2A–D. Cell wall protein selection procedure and statistical data. Supplementary Figure S1: Representative gel image for each of the fractions (CaCl2, EGTA, and LiCl). Supplementary Figure S2: Sequence alignment of the M. truncatula NCBInr gi: 922395483, the PNGase identified in this study with its homologues NCBInr gi:56405352 from almonds, and NCBInr gi:723633692 from tomatoes. Alignment was carried out with Clustal Omega http://www.ebi.ac.uk/Tools/msa/clustalo/. Data Availability: The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (https://www.ebi.ac.uk/pride/archive/) with the dataset identifier PXD010666.

Author Contributions

B.P., K.S., S.L., and J.-F.H. conceived the study and the experiment. B.P., K.S., G.G., and J.R. performed the analyses and did the interpretation of the data. B.P., G.G., and K.S. wrote the manuscript which was read and edited by all co-authors that approved the manuscript for submission.

Funding

This work was supported by the National Research Fund, Luxembourg [FNR Project CADWALL (FNR/FWO project INTER/FWO/12/14)].

Acknowledgments

We acknowledge Céline Leclercq, Sébastien Planchon, and Laurent Solinhac for their day-to-day technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guerriero, G.; Hausman, J.F.; Cai, G. No Stress! Relax! Mechanisms governing growth and shape in plant cells. Int. J. Mol. Sci. 2014, 15, 5094–5114. [Google Scholar] [CrossRef] [PubMed]
  2. Underwood, W. The plant cell wall: A dynamic barrier against pathogen invasion. Front. Plant Sci. 2012, 3, 85. [Google Scholar] [CrossRef] [PubMed]
  3. Cosgrove, D.J. Plant cell wall extensibility: Connecting plant cell growth with cell wall structure, mechanics, and the action of wall-modifying enzymes. J. Exp. Bot. 2016, 67, 463–476. [Google Scholar] [CrossRef] [PubMed]
  4. Houston, K.; Tucker, M.R.; Chowdhury, J.; Shirley, N.; Little, A. The plant cell wall: A complex and dynamic structure as revealed by the responses of genes under stress conditions. Front. Plant Sci. 2016, 7, 984. [Google Scholar] [CrossRef] [PubMed]
  5. Cassab, G.I.; Varner, J.E. Cell wall proteins. Annu. Rev. Plant Physiol. Plant Mol. Biol. 1988, 39, 321–353. [Google Scholar] [CrossRef]
  6. Boaretto, L.F.; Mazzafera, P. The proteomes of feedstocks used for the production of second-generation ethanol: A lacuna in the biofuel era. Ann. Appl. Biol. 2013, 163, 12–22. [Google Scholar] [CrossRef]
  7. Rose, J.K.C.; Lee, S.J. Straying off the highway: Trafficking of secreted plant proteins and complexity in the plant cell wall proteome. Plant Physiol. 2010, 153, 433–436. [Google Scholar] [CrossRef]
  8. Pinedo, M.; Regente, M.; Elizande, M.; Quiroga, I.Y.; Pagnussat, A.L.; Jorrin-Novo, J.; Maldonado, A.; de la Canal, L. Extracellular sunflower proteins: Evidence on non-classical secretion of a jacalin-related lectin. Protein Pep. Lett. 2012, 19, 270–276. [Google Scholar] [CrossRef]
  9. Rayon, C.; Lerouge, P.; Faye, L. The protein N-glycosylation in plants. J. Exp. Bot. 1998, 49, 1463–1472. [Google Scholar] [CrossRef]
  10. Hijazi, M.; Velasquez, S.M.; Jamet, E.; Estevez, J.M.; Albenne, C. An update on post-translational modifications of hydroxyproline-rich glycoproteins: Toward a model highlighting their contribution to plant cell wall architecture. Front. Plant Sci. 2014, 5, 395. [Google Scholar] [CrossRef]
  11. de la Canal, L.; Pinedo, M. Extracellular vesicles: A missing component in plant cell wall remodeling. J. Exp. Bot. 2018, 69, 4655–4658. [Google Scholar] [CrossRef] [PubMed]
  12. Yun, H.S.; Kwon, C.; Kang, B.G.; Lee, J.S.; Han, T.J.; Chang, S.C.; Kim, S.K. A xyloglucan endotransglycosylase/hydrolase1, VrXTH1, is associated with cell elongation in mungbean hypocotyls. Physiol. Plant. 2005, 125, 106–113. [Google Scholar] [CrossRef]
  13. Micheli, F. Pectin methylesterases: Cell wall enzymes with important roles in plant physiology. Trends Plant Sci. 2001, 6, 414–419. [Google Scholar] [CrossRef]
  14. Altmann, F.; Paschinger, K.; Dalik, T.; Vorauer, K. Characterisation of peptide-N4-(N-acetyl-β-glucosaminyl)asparagine amidase A and its N-glycans. Eur. J. Biochem. 1998, 252, 118–123. [Google Scholar] [CrossRef] [PubMed]
  15. Gantulga, D.; Ahn, Y.O.; Zhou, C.; Battogtokh, D.; Bevan, D.R.; Winkel, B.S.J.; Esen, A. Comparative characterization of the Arabidopsis subfamily a1 β-galactosidases. Phytochemistry 2009, 70, 1999–2009. [Google Scholar] [CrossRef]
  16. Vazquez-Lobo, A.; Roujol, D.; Zuniga-Sanchez, E.; Albenne, C.; Pinero, D.; Gamboa deBuen, A.; Jamet, E. The highly conserved spermatophyte cell wall DUF642 protein family: Phylogeny and first evidence of interaction with cell wall polysaccharides in vitro. Mol. Phylogenet. Evol. 2012, 63, 510–520. [Google Scholar] [CrossRef]
  17. Zhao, Q.; Nakashima, J.; Chen, F.; Yin, Y.; Fu, C.; Yun, J.; Shao, H.; Wang, X.; Wang, Z.Y.; Dixon, R.A. LACCASE is necessary and nonredundant with PEROXIDASE for lignin polymerization during vascular development in Arabidopsis. Plant Cell 2013, 25, 3976–3987. [Google Scholar] [CrossRef]
  18. Zúñiga-Sánchez, E.; Gamboa-deBuen, A. The two DUF642 At5g11420 and At4g32460-encoded proteins interact in vitro with the AtPME3 catalytic domain. In Protein Interactions; Cai, J., Wang, R.E., Eds.; InTechOpen: Rijeka, Croatia, 2012; Chapter 7; p. 464. [Google Scholar] [CrossRef]
  19. Verdonk, J.C.; Hatfield, R.D.; Sullivan, M.L. Proteomic analysis of cell walls of two developmental stages of alfalfa stems. Front. Plant Sci. 2012, 3, 279. [Google Scholar] [CrossRef]
  20. Feiz, L.; Irshad, M.; Pont-Lezica, R.F.; Canut, H.; Jamet, E. Evaluation of cell wall preparations for proteomics: A new procedure for purifying cell walls from Arabidopsis hypocotyls. Plant Methods 2006, 2, 10. [Google Scholar] [CrossRef]
  21. Printz, B.; Dos Santos Morais, R.; Wienkoop, S.; Sergeant, K.; Lutts, S.; Hausman, J.F.; Renaut, J. An improved protocol to study the plant cell wall proteome. Front. Plant Sci. 2015, 6, 237. [Google Scholar] [CrossRef] [Green Version]
  22. Barros, J.; Temple, S.; Dixon, R.A. Development and commercialization of reduced lignin alfalfa. Curr. Opin. Biotechnol. 2019, 56, 48–54. [Google Scholar] [CrossRef]
  23. Gutsch, A.; Keunen, E.; Guerriero, G.; Renaut, J.; Cuypers, A.; Hausman, J.F.; Sergeant, K. Long-term cadmium exposure influences the abundance of proteins that impact the cell wall structure in Medicago sativa stems. Plant Biol. 2018, 20, 1023–1035. [Google Scholar] [CrossRef]
  24. Gutsch, A.; Sergeant, K.; Keunen, E.; Prinsen, E.; Guerriero, G.; Renaut, J.; Hausman, J.F.; Cuypers, A. Does long-term cadmium exposure influence the composition of pectic polysaccharides in the cell wall of Medicago sativa stems? BMC Plant Biol. 2019, 19, 271. [Google Scholar] [CrossRef]
  25. Printz, B.; Guerriero, G.; Sergeant, K.; Renaut, J.; Lutts, S.; Hausman, J.F. Ups and downs in alfalfa: Proteomic and metabolic changes occurring in the growing stem. Plant Sci. 2015, 238, 13–25. [Google Scholar] [CrossRef] [Green Version]
  26. Duruflé, H.; San Clemente, H.; Balliau, T.; Zivy, M.; Dunand, C.; Jamet, E. Cell wall proteome analysis of Arabidopsis thaliana mature stems. Proteomics 2017, 17, 1600449. [Google Scholar] [CrossRef]
  27. Calderan-Rodrigues, M.J.; Jamet, E.; Douché, T.; Bonassi, M.B.; Cataldi, T.R.; Fonseca, J.G.; San Clemente, H.; Pont-Lezica, R.; Labate, C.A. Cell wall proteome of sugarcane stems: Comparison of a destructive and a non-destructive extraction method showed differences in glycoside hydrolases and peroxidases. BMC Plant Biol. 2016, 16, 14. [Google Scholar] [CrossRef]
  28. Chabi, M.; Goulas, E.; Leclercq, C.C.; de Waele, I.; Rihouey, C.; Cenci, U.; Day, A.; Blervacq, A.S.; Neutelings, G.; Duponchel, L.; et al. A cell wall proteome and targeted cell wall analyses provide novel information on hemicellulose metabolism in flax. Mol. Cell. Proteom. 2017, 16, 1634–1651. [Google Scholar] [CrossRef]
  29. Komatsu, S.; Yanagawa, Y. Cell wall proteomics of crops. Front. Plant Sci. 2013, 4, 17. [Google Scholar] [CrossRef] [Green Version]
  30. Sergeant, K.; Printz, B.; Gutsch, A.; Behr, M.; Renaut, J.; Hausman, J.F. Didehydrophenylalanine, an abundant modification in the beta subunit of plant polygalacturonases. PLoS ONE 2017, 12, e0171990. [Google Scholar] [CrossRef]
  31. Shpak, E.; Leykam, J.F.; Kieliszewski, M.J. Synthetic genes for glycoprotein design and the elucidation of hydroxyproline-O-glycosylation codes. Proc. Natl. Acad. Sci. USA 1999, 96, 14736–14741. [Google Scholar] [CrossRef]
  32. Hoagland, D.R.; Arnon, D.I. The water-culture method for growing plants without soil. Circ. Calif. Agric. Exp. Stn. 1950, 347, 32. [Google Scholar]
  33. Printz, B.; Sergeant, K.; Lutts, S.; Guignard, C.; Renaut, J.; Hausman, J.F. From tolerance to acute metabolic deregulation: Contribution of proteomics to dig into the molecular response of alder-species under a polymetallic exposure. J. Proteome Res. 2013, 12, 5160–5179. [Google Scholar] [CrossRef]
  34. Almagro Armenteros, J.J.; Sønderby, C.K.; Sønderby, S.K.; Nielsen, H.; Winther, O. DeepLoc: Prediction of protein subcellular localization using deep learning. Bioinformatics 2017, 33, 3387–3395. [Google Scholar] [CrossRef]
  35. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R package for multivariate analysis. J. Stat. Softw. 2008, 25, 1–18. [Google Scholar] [CrossRef]
  36. Fawal, N.; Li, Q.; Savelli, B.; Brette, M.; Passaia, G.; Fabre, M.; Mathé, C.; Dunand, C. PeroxiBase: A database for large-scale evolutionary analysis of peroxidases. Nucleic Acids Res. 2013, 41, D441–D444. [Google Scholar] [CrossRef]
  37. Watson, B.S.; Lei, Z.; Dixon, R.A.; Sumner, L.W. Proteomics of Medicago sativa cell walls. Phytochemistry 2004, 65, 1709–1720. [Google Scholar] [CrossRef]
  38. Irshad, M.; Canut, H.; Borderies, G.; Pont-Lezica, R.; Jamet, E. A new picture of cell wall protein dynamics in elongating cells of Arabidopsis thaliana: Confirmed actors and newcomers. BMC Plant Biol. 2008, 8, 94. [Google Scholar] [CrossRef]
  39. Douché, T.; San Clemente, H.; Burlat, V.; Roujol, D.; Valot, B.; Zivy, M.; Pont-Lezica, R.; Jamet, E. Brachypodium distachyon as a model plant toward improved biofuel crops: Search for secreted proteins involved in biogenesis and disassembly of cell wall polymers. Proteomics 2013, 13, 2438–2454. [Google Scholar] [CrossRef]
  40. Emanuelsson, O.; Brunak, S.; von Heijne, G.; Nielsen, H. Locating proteins in the cell using TargetP, SignalP and related tools. Nat. Protoc. 2007, 2, 953–971. [Google Scholar] [CrossRef]
  41. Neuhaus, J.M.; Pietrzak, M.; Boller, T. Mutation analysis of the C-terminal vacuolar targeting peptide of tobacco chitinase: Low specificity of the sorting system, and gradual transition between intracellular retention and secretion into the extracellular space. Plant J. 1994, 5, 45–54. [Google Scholar] [CrossRef]
  42. Jamet, E.; Albenne, C.; Boudart, G.; Irshad, M.; Canut, H.; Pont-Lezica, R. Recent advances in plant cell wall proteomics. Proteomics 2008, 8, 893–908. [Google Scholar] [CrossRef] [Green Version]
  43. Wasąg, P.; Suwińska, A.; Zakrzewski, P.; Walczewski, J.; Lenartowski, R.; Lenartowska, M. Calreticulin localizes to plant intra/extracellular peripheries of highly specialized cells involved in pollen-pistil interactions. Protoplasma 2018, 255, 57–67. [Google Scholar] [CrossRef]
  44. Jia, X.Y.; He, L.H.; Jing, R.L.; Li, R.Z. Calreticulin: Conserved protein and diverse functions in plants. Physiol. Plant. 2009, 136, 127–138. [Google Scholar] [CrossRef]
  45. Hojrup, P.; Roepstorff, P.; Houen, G. Human placental calreticulin characterization of domain structure and post-translational modifications. Eur. J. Biochem. 2001, 268, 2558–2565. [Google Scholar] [CrossRef]
  46. Cai, Y.-M.; Yu, J.; Ge, Y.; Mironov, A.; Gallois, P. Two proteases with caspase-3-like activity, cathepsin B and proteasome, antagonistically control ER-stress-induced programmed cell death in Arabidopsis. New Phytol. 2018, 218, 1143–1155. [Google Scholar] [CrossRef]
  47. Tarr, M.; Young, P.J.; Morse, R.; Shaw, D.J.; Haigh, R.; Petrov, P.G.; Johnson, S.J.; Winyard, P.G.; Eggleton, P. A mechanism of release of calreticulin from cells during apoptosis. J. Mol. Biol. 2010, 401, 799–812. [Google Scholar] [CrossRef]
  48. Toby, T.K.; Fornelli, L.; Kelleher, N.L. Progress in top-down proteomics and the analysis of proteoforms. Annu. Rev. Anal. Chem. 2016, 9, 499–519. [Google Scholar] [CrossRef]
  49. Printz, B.; Guerriero, G.; Sergeant, K.; Audinot, J.-N.; Guignard, C.; Renaut, J.; Lutts, S.; Hausman, J.F. Combining-omics to unravel the impact of copper nutrition on alfalfa (Medicago sativa) stem metabolism. Plant Cell Physiol. 2016, 57, 407–422. [Google Scholar] [CrossRef]
  50. O’Rourke, J.A.; Fu, F.; Bucciarelli, B.; Yang, S.S.; Samac, D.A.; Lamb, J.F.S.; Monteros, M.J.; Graham, M.A.; Gronwald, J.W.; Krom, N.; et al. The Medicago sativa gene index 1.2: A web-accessible gene expression atlas for investigating expression differences between Medicago sativa subspecies. BMC Genom. 2015, 16, 502. [Google Scholar] [CrossRef]
  51. Watson, C.F.; Zheng, L.; DellaPenna, D. Reduction of tomato polygalacturonase beta subunit expression affects pectin solubilization and degradation during fruit ripening. Plant Cell 1994, 6, 1623–1634. [Google Scholar] [CrossRef]
  52. Orozco-Cárdenas, M.L.; Ryan, C.A. Polygalacturonase β-Subunit antisense gene expression in tomato plants leads to a progressive enhanced wound response and necrosis in leaves and abscission of developing flowers. Plant Physiol. 2003, 133, 693–701. [Google Scholar] [CrossRef]
  53. Yuan, Z.; Yao, X.; Zhang, D.; Sun, Y.; Huang, H. Genome-wide expression profiling in seedlings of the Arabidopsis mutant uro that is defective in the secondary cell wall formation. J. Integr. Plant Biol. 2007, 49, 1754–1762. [Google Scholar] [CrossRef]
  54. Yun, H.S.; Kwon, C.; Kim, T.W.; Joo, S.H.; Cho, M.H.; Kang, B.G.; Chang, S.C.; Kim, S.K. Regulation of VrXTH1 expression in mungbean. J. Plant Biol. 2007, 50, 65–69. [Google Scholar] [CrossRef]
  55. Campbell, P.; Braam, J. In vitro activities of four xyloglucan endotransglycosylases from Arabidopsis. Plant J. 2002, 18, 371–382. [Google Scholar] [CrossRef]
  56. Morales-Quintana, L.; Carrasco-Orellana, C.; Beltran, D.; Moya-Leon, M.A.; Herrera, R. Molecular insights of a xyloglucan endo-transglycosylase/hydrolase of radiata pine (PrXTH1) expressed in response to inclination: Kinetics and computational study. Plant Physiol. Biochem. 2019, 136, 155–161. [Google Scholar] [CrossRef]
  57. Salazar-Iribe, A.; Cruz-Valderrama, J.E.; Jimenez-Duran, K.; Gomez-Maqueo, X.; Gamboa-deBuen, A. BIDDXI, a DUF642 cell wall protein, is involved in hypocotyl growth via auxin efflux. J. Plant Physiol. 2018, 231, 105–109. [Google Scholar] [CrossRef]
  58. Zúñiga-Sánchez, E.; Soriano, D.; Martínez-Barajas, E.; Orozco-Segovia, A.; Gamboa-deBuen, A. BIIDXI, the At4g32460 DUF642 gene, is involved in pectin methyl esterase regulation during Arabidopsis thaliana seed germination and plant development. BMC Plant Biol. 2014, 14, 338. [Google Scholar] [CrossRef]
  59. Dekkers, B.J.W.; Pearce, S.; van Bolderen-Veldkamp, R.P.; Marshall, A.; Widera, P.; Gilbert, J.; Drost, H.G.; Bassel, G.W.; Müller, K.; King, J.R.; et al. Transcriptional dynamics of two seed compartments with opposing roles in Arabidopsis seed germination. Plant Physiol. 2013, 163, 205–215. [Google Scholar] [CrossRef]
  60. Araki, R.; Mermod, M.; Yamasaki, H.; Kamiya, T.; Fuijwara, T.; Shikanai, T. SPL7 locally regulates copper-homeostasis-related genes in Arabidopsis. J. Plant Physiol. 2018, 224–225, 137–143. [Google Scholar] [CrossRef]
  61. Xu, P.; Cai, X.T.; Wang, Y.; Xing, L.; Chen, Q.; Xiang, C.-B. HDG11 upregulates cell-wall-loosening protein genes to promote root elongation in Arabidopsis. J. Exp. Bot. 2014, 65, 4285–4295. [Google Scholar] [CrossRef]
  62. Jung, J.-H.; Ju, Y.; Seo, P.J.; Lee, J.-H.; Park, C.-M. The SOC1-SPL module integrates photoperiod and gibberellic acid signals to control flowering time in Arabidopsis. Plant J. Cell. Mol. Biol. 2012, 69, 577–588. [Google Scholar] [CrossRef]
  63. Salazar-Iribe, A.; Agredano-Moreno, L.T.; Zúñiga-Sánchez, E.; Jiménez-Garcia, L.F.; Gamboa-deBuen, A. The cell wall DUF642 At2g41800 (TEB) protein is involved in hypocotyl cell elongation. Plant Sci. 2016, 253, 206–214. [Google Scholar] [CrossRef]
  64. Iglesias, N.; Abelenda, J.A.; Rodiño, M.; Sampedro, J.; Revilla, G.; Zarra, I. Apoplastic glycosidases active against xyloglucan oligosaccharides of Arabidopsis thaliana. Plant Cell Physiol. 2006, 47, 55–63. [Google Scholar] [CrossRef]
  65. Tesfaye, M.; Yang, S.S.; Lamb, J.F.S.; Jung, H.J.G.; Samac, D.A.; Vance, C.P.; Gronwald, J.W.; VandenBosch, K.A. Medicago truncatula as a model for dicot cell wall development. BioEnergy Res. 2009, 2, 59–76. [Google Scholar] [CrossRef]
  66. Cosio, C.; Dunand, C. Specific functions of individual class III peroxidase genes. J. Exp. Bot. 2009, 60, 391–408. [Google Scholar] [CrossRef]
  67. Francoz, E.; Ranocha, P.; Nguyen-Kim, H.; Jamet, E.; Burlat, V.; Dunand, C. Roles of cell wall peroxidases in plant development. Phytochemistry 2014, 112, 15–21. [Google Scholar] [CrossRef]
  68. Behr, M.; Legay, S.; Hausman, J.-F.; Guerriero, G. Analysis of cell wall-related genes in organs of Medicago sativa L. under different abiotic stresses. Int. J. Mol. Sci. 2015, 16, 16104–16124. [Google Scholar] [CrossRef]
  69. Rahman, M.A.; Alam, I.; Kim, Y.-G.; Ahn, N.Y.; Heo, S.H.; Lee, D.G.; Liu, G.; Lee, B.H. Screening for salt-responsive proteins in two contrasting alfalfa cultivars using a comparative proteome approach. Plant Physiol. Biochem. 2015, 89, 112–122. [Google Scholar] [CrossRef]
  70. Soares, N.C.; Wojtkowska, J.; Jackson, P.A. A proteomic analysis of the wound response in Medicago leaves reveals the early activation of a ROS-sensitive signal pathway. J. Proteom. 2011, 74, 1411–1420. [Google Scholar] [CrossRef]
  71. Oger, E.; Marino, D.; Guigonis, J.-M.; Pauly, N.; Puppo, A. Sulfenylated proteins in the Medicago truncatula–Sinorhizobium meliloti symbiosis. J. Proteom. 2012, 75, 4102–4113. [Google Scholar] [CrossRef]
  72. Li, N.; Zhao, M.; Liu, T.; Dong, L.; Cheng, Q.; Wu, J.; Wang, L.; Chen, X.; Zhang, C.; Lu, W.; et al. A novel soybean dirigent gene GmDIR22 contributes to promotion of lignan biosynthesis and enhances resistance to Phytophthora sojae. Front. Plant Sci. 2017, 8, 1185. [Google Scholar] [CrossRef]
  73. Bacete, L.; Melida, H.; Miedes, E.; Molina, A. Plant cell wall-mediated immunity: Cell wall changes trigger disease resistance responses. Plant J. 2018, 93, 614–636. [Google Scholar] [CrossRef]
  74. Song, M.; Peng, X. Genome-wide identification and characterization of DIR genes in Medicago truncatula. Biochem. Genet. 2019, 57, 1–20. [Google Scholar] [CrossRef]
  75. Francin-Allami, M.; Lollier, V.; Pavlovic, M.; San Clemente, H.; Rogniaux, H.; Jamet, E.; Guillon, F.; Larré, C. Understanding the remodelling of cell walls during Brachypodium distachyon grain development through a sub-cellular quantitative proteomic approach. Proteomes 2016, 4, 21. [Google Scholar] [CrossRef]
  76. Tarentino, A.L.; Gomez, C.M.; Plummer, T.H. Deglycosylation of asparagine-linked glycans by peptide:N-glycosidase F. Biochemistry 1985, 24, 4665–4671. [Google Scholar] [CrossRef]
  77. Maeda, M.; Kimura, Y. Structural features of free N-glycans occurring in plants and functional features of de-N-glycosylation enzymes, ENGase, and PNGase: The presence of unusual plant complex type N-glycans. Front. Plant Sci. 2014, 5, 429. [Google Scholar] [CrossRef]
  78. Boudart, G.; Jamet, E.; Rossignol, M.; Lafitte, C.; Borderies, G.; Jauneau, A.; Esquerré-Tugayé, M.T.; Pont-Lezica, R. Cell wall proteins in apoplastic fluids of Arabidopsis thaliana rosettes: Identification by mass spectrometry and bioinformatics. Proteomics 2005, 5, 212–221. [Google Scholar] [CrossRef]
  79. Tsai, A.Y.; Kunieda, T.; Rogalski, J.; Foster, L.J.; Ellis, B.E.; Haughn, G.W. Identification and characterization of Arabidopsis seed coat mucilage proteins. Plant Physiol. 2017, 173, 1059–1074. [Google Scholar] [CrossRef]
  80. Charmont, S.; Jamet, E.; Pont-Leziva, R.; Canut, H. Proteomics analysis of secreted proteins from Arabidopsis thaliana seedlings: Improved recovery following removal of phenolic compounds. Phytochemistry 2004, 66, 453–461. [Google Scholar] [CrossRef]
  81. Gu, C.; Shabab, M.; Strasser, R.; Wolters, P.J.; Shindo, T.; Niemer, M.; Kaschani, F.; Mach, L.; van der Hoorn, R.A. Post-translational regulation and trafficking of the granulin-containing protease RD21 of Arabidopsis thaliana. PLoS ONE 2012, 7, e32422. [Google Scholar] [CrossRef]
  82. Jashni, M.K.; Mehrabi, R.; Collemare, J.; Mesarich, C.H.; de Wit, P.J.G.M. The battle in the apoplast: Further insights into the roles of proteases and their inhibitors in plant-pathogen interactions. Front. Plant Sci. 2015, 6, 584. [Google Scholar] [CrossRef]
  83. Pogorelko, G.V.; Juvale, P.S.; Rutter, W.B.; Hütten, M.; Maier, T.R.; Hewezi, T.; Paulus, J.; van der Hoorn, R.A.L.; Grundler, F.M.W.; Siddique, S.; et al. Re-targeting of a plant defense protease by a cyst nematode effector. Plant J. 2019, 98, 1000–1014. [Google Scholar] [CrossRef]
  84. Cvitanich, C.; Przybylowicz, W.J.; Urbanski, D.F.; Jurkiewicz, A.M.; Mesjasz-Przybylowicz, J.; Blair, M.W.; Astudillo, C.; Jensen, E.Ø.; Stougaard, J. Iron and ferritin accumulate in separate cellular locations in Phaseolus seeds. BMC Plant Biol. 2010, 10, 26. [Google Scholar] [CrossRef]
  85. Ravet, K.; Touraine, B.; Boucherez, J.; Briat, J.F.; Gaymard, F.; Cellier, F. Ferritins control interaction between iron homeostasis and oxidative stress in Arabidopsis. Plant J. 2009, 57, 400–412. [Google Scholar] [CrossRef]
  86. Strozycki, P.M.; Szczurek, A.; Lotocka, B.; Figlerowicz, M.; Legocki, A.B. Ferritins and nodulation in Lupinus luteus: Iron management in indeterminate type nodules. J. Exp. Bot. 2007, 58, 3145–3153. [Google Scholar] [CrossRef]
  87. Hotte, N.; Deyholos, M. A flax fibre proteome: Identification of proteins enriched in bast fibres. BMC Plant Biol. 2008, 8, 52. [Google Scholar] [CrossRef]
  88. De Hoff, P.L.; Brill, L.M.; Hirsch, A.M. Plant lectins: The ties that bind in root symbiosis and plant defense. Mol. Genet. Genom. 2009, 282, 1–15. [Google Scholar] [CrossRef]
  89. Dias, R.O.; Machado, L.S.; Migliolo, L.; Franco, O.L. Insights into animal and plant lectins with antimicrobial activities. Molecules 2015, 20, 519–541. [Google Scholar] [CrossRef]
Figure 1. Principal component analysis (PCA) showing the distribution of the individuals based on the set of 113 proteins differentially regulated in the different stem regions. Only the two first components are represented (PC1: 55.98%, PC2: 10.24%). Each dot refers to one replicate and squares represent the barycenter of the dots corresponding to one stem region.
Figure 1. Principal component analysis (PCA) showing the distribution of the individuals based on the set of 113 proteins differentially regulated in the different stem regions. Only the two first components are represented (PC1: 55.98%, PC2: 10.24%). Each dot refers to one replicate and squares represent the barycenter of the dots corresponding to one stem region.
Biology 08 00060 g001
Figure 2. Profiling of the CWP and non-CWP (intracellular) discussed in this study and identified in the apical, intermediate, and basal stem region (big font size = high abundance, medium font size = medium abundance).
Figure 2. Profiling of the CWP and non-CWP (intracellular) discussed in this study and identified in the apical, intermediate, and basal stem region (big font size = high abundance, medium font size = medium abundance).
Biology 08 00060 g002
Figure 3. Conserved motifs in the promoters of the seven DUF642 genes from M. truncatula. The 1000 bp upstream of the start codon were analyzed using the MEME suite 4.11.2. The first motif in each table cell is the one found in the database; the second motif is the one found in the promoters of alfalfa genes. The p-values are indicated.
Figure 3. Conserved motifs in the promoters of the seven DUF642 genes from M. truncatula. The 1000 bp upstream of the start codon were analyzed using the MEME suite 4.11.2. The first motif in each table cell is the one found in the database; the second motif is the one found in the promoters of alfalfa genes. The p-values are indicated.
Biology 08 00060 g003
Table 1. Clustering of the proteins according to the procedure described in the method section: By stem region (basal part (Base); intermediate part (Int); apical part (Apex) of the stem), the average abundance measured for each spot is given as log(2) of the standardized volume and represented by a green bar. The column “Loc” and “Type” respectively indicate the localization and protein type (soluble (Sol) and membranair (Memb)), as predicted by Deeploc. Extracellular (Extrac); cytoplasmic (Cytop); plastidial (Plast); Lysosomal/vacuolar (Lyso). Detailed identification data are given in Supplementary Table S1A for the CaCl2 fraction, Supplementary Table S1B for the EGTA fraction, and Supplementary Table S1C for the LiCl fraction. The functional classification of the spots that contribute significantly to two PCA dimensions are given in Supplementary Table S2D.
Table 1. Clustering of the proteins according to the procedure described in the method section: By stem region (basal part (Base); intermediate part (Int); apical part (Apex) of the stem), the average abundance measured for each spot is given as log(2) of the standardized volume and represented by a green bar. The column “Loc” and “Type” respectively indicate the localization and protein type (soluble (Sol) and membranair (Memb)), as predicted by Deeploc. Extracellular (Extrac); cytoplasmic (Cytop); plastidial (Plast); Lysosomal/vacuolar (Lyso). Detailed identification data are given in Supplementary Table S1A for the CaCl2 fraction, Supplementary Table S1B for the EGTA fraction, and Supplementary Table S1C for the LiCl fraction. The functional classification of the spots that contribute significantly to two PCA dimensions are given in Supplementary Table S2D.
ClusterSpotgi-NumberNameLocTypeBaseInt.Apex
1CaCl2-2127gi|2951684lectin [M. sativa]ExtracSol0.37−0.20−0.39
1CaCl2-2128gi|2951684lectin [M. sativa]ExtracSol0.25−0.30−0.27
1CaCl2-2193gi|400180RecName: Full=Truncated lectin 2ExtracSol0.34−0.14−0.14
1EGTA-1167gi|388505450Glutelin type-A [M. truncatula]CytopSol0.25−0.14−0.23
1EGTA-1390gi|400180RecName: Full=Truncated lectin 2ExtracSol0.40−0.20−0.23
1EGTA-1411gi|2951684lectin [M. sativa]ExtracSol0.20−0.29−0.30
1EGTA-1420gi|400180RecName: Full=Truncated lectin 2ExtracSol0.19−0.19−0.22
1LiCl-765gi|2951684lectin [M. sativa]ExtracSol0.21−0.25−0.23
2CaCl2-1870gi|388505450Glutelin type-A [M. truncatula]CytopSol0.20−0.20−0.26
2CaCl2-2363gi|357468557Ferritin-3 [M. truncatula]PlastSol0.24−0.04−0.44
2CaCl2-2627gi|357492793Ferritin-2 [M. truncatula]PlastSol0.11−0.01−0.40
2EGTA-2050gi|50320305putative defensin 2.1 precursor [M. sativa]ExtracSol0.15−0.14−0.19
3CaCl2-1013gi|922395483PNGase A [M. truncatula]ExtracSol0.02−0.08−0.17
3CaCl2-2165gi|228204925chitinase class III-1 [M. sativa]ExtracSol0.04−0.12−0.21
3CaCl2-2426gi|357468557Ferritin-3 [M. truncatula]PlastSol0.09−0.08−0.14
3EGTA-1228gi|357437719Cysteine proteinase [M. truncatula]LysoSol0.16−0.12−0.14
3EGTA-1372gi|228204925chitinase class III-1 [M. sativa]ExtracSol0.03−0.14−0.23
3EGTA-2117gi|3913031RecName: Full=Beta-amylase [M. sativa]CytopSol0.12−0.11−0.20
3EGTA-570gi|3913031RecName: Full=Beta-amylase [M. sativa]CytopSol0.11−0.10−0.22
3EGTA-573gi|3913031RecName: Full=Beta-amylase [M. sativa]CytopSol0.09−0.06−0.19
3EGTA-586gi|3913031RecName: Full=Beta-amylase [M. sativa]CytopSol0.12−0.09−0.18
3EGTA-587gi|3913031RecName: Full=Beta-amylase [M. sativa]CytopSol0.07−0.05−0.14
3EGTA-593gi|3913031RecName: Full=Beta-amylase [M. sativa]CytopSol0.08−0.13−0.22
3EGTA-595gi|3913031RecName: Full=Beta-amylase [M. sativa]CytopSol0.05−0.11−0.16
3EGTA-610gi|3913031RecName: Full=Beta-amylase [M. sativa]CytopSol0.12−0.09−0.28
3EGTA-615gi|3913031RecName: Full=Beta-amylase [M. sativa]CytopSol0.10−0.15−0.23
3EGTA-640gi|3913032RecName: Full=Beta-amylase [M. sativa]CytopSol0.13−0.07−0.14
4CaCl2-1716gi|388518933Endochitinase [M. truncatula]ExtracSol0.260.10−0.12
4CaCl2-1880gi|357448991endo-1,3-beta-glucosidase [M. truncatula]ExtracSol0.190.06−0.08
4CaCl2-2356gi|357468557Ferritin-3 [M. truncatula]PlastSol0.140.08−0.20
4CaCl2-2425gi|357468557Ferritin-3 [M. truncatula]PlastSol0.230.06−0.32
4CaCl2-2439gi|357468557Ferritin-3 [M. truncatula]PlastSol0.280.04−0.21
4CaCl2-2553gi|357468557Ferritin-3 [M. truncatula]PlastSol0.24−0.04−0.18
4CaCl2-2827gi|22901738high FW vegetative storage proteinExtracSol0.21−0.03−0.23
4CaCl2-2828gi|357437721Cysteine proteinase [M. truncatula]LysoSol0.090.05−0.13
4EGTA-1111gi|71534882serine carboxypeptidase S10 [M. sativa]ExtracSol0.210.00−0.16
4EGTA-1114gi|3913031RecName: Full=Beta-amylase [M. sativa]CytopSol0.180.06−0.08
4EGTA-1119gi|50317234Aspartic proteinase Asp1 [M. truncatula]ExtracSol0.170.00−0.12
4EGTA-1143gi|357448997endo-1,3-beta-glucosidase [M. truncatula]ExtracSol0.15−0.02−0.10
4EGTA-1166gi|298364452chitinase [M. sativa]ExtracSol0.110.03−0.17
4EGTA-1168gi|657377089peroxidase family protein [M. truncatula]ExtracSol0.16−0.03−0.12
4EGTA-1233gi|357437719Cysteine proteinase [M. truncatula]ExtracSol0.120.03−0.14
4EGTA-1461gi|169147017putative thaumatin-like proteinExtracSol0.190.05−0.04
4EGTA-1922gi|71534922thioredoxin h [M. sativa]CytopSol0.12−0.01−0.15
4EGTA-914gi|657371151eukaryotic aspartyl proteaseExtracSol0.19−0.02−0.17
4LiCl-1068gi|23049900put early nodulin-like 2 [T pratense]CMMemb0.100.09−0.08
5CaCl2-1769gi|357491415Peroxidase [M. truncatula]ExtracSol0.400.290.11
5CaCl2-2269gi|3688528pre-pro-TPE4A protein [Pisum sativum]LysoSol−0.010.170.05
5EGTA-1052gi|357517805Endochitinase [M. truncatula]LysoSol0.180.070.01
5EGTA-1057gi|388494834Endochitinase [M. truncatula]LysoSol0.230.09−0.01
5EGTA-1450gi|3688528cysteine proteinase [Vicia sativa]LysoSol0.000.210.06
5EGTA-1451gi|571536058PREDICTED: vignain-like [Glycine max]LysoSol−0.020.150.00
5EGTA-1452gi|30141021cysteine protease-2 [Helianthus annuus]LysoSol−0.040.14−0.01
5EGTA-2122gi|357448991endo-1,3-beta-glucosidase [M. truncatula]ExtracSol0.260.10−0.03
6LiCl-885gi|357513969Germin-like protein [M. truncatula]ExtracSol−0.40−0.28−0.08
6LiCl-935gi|388502800Germin-like protein [Pisum sativum]ExtracSol−0.31−0.25−0.14
7CaCl2-1593gi|593701343Cyclophilin 38 isoform 1 [T. cacao]PlastMemb−0.15−0.060.13
7EGTA-1388gi|357480321Lactoylglutathione lyase [M. truncatula]CytopSol−0.15−0.020.04
7EGTA-1480gi|357512271Ribose-5-phosphate isomer. [M. truncatula]PlastSol−0.180.010.18
7EGTA-1586gi|388502800Germin-like protein [Pisum sativum]ExtracSol−0.23−0.07−0.01
7EGTA-196gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.180.040.01
7EGTA-197gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.190.030.03
7EGTA-201gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.240.020.03
7EGTA-203gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.230.050.05
7EGTA-205gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.190.010.02
7EGTA-206gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.190.020.03
7EGTA-2120gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.160.050.03
7EGTA-214gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.240.040.03
7EGTA-215gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.230.050.05
7EGTA-216gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.180.020.01
7EGTA-220gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.210.060.06
7EGTA-221gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.180.060.06
7EGTA-222gi|657400954beta-like galactosidase [M. truncatula]ExtracSol−0.230.000.04
7EGTA-240gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.230.020.01
7EGTA-272gi|87241408Protease-associated PA [M. truncatula]ExtracSol−0.08−0.010.14
7EGTA-326gi|357437609Beta-galactosidase [M. truncatula]ExtracSol-0.220.040.02
7EGTA-333gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.160.020.06
7EGTA-341gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.160.030.04
7EGTA-342gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.190.050.01
7EGTA-379gi|357437609Beta-galactosidase [M. truncatula]ExtracSol−0.20−0.02−0.05
8CaCl2-2612gi|357476945Acid beta-fructofuranosidase LysoSol−0.160.020.27
8CaCl2-2968gi|388514479Oxygen-evolving enhancer protein 2PlastSol−0.16-0.030.28
8EGTA-1296gi|357438145Cysteine proteinase [M. truncatula]LysoSol−0.040.030.21
8EGTA-1485gi|357512271Ribose-5-phosphate isomer. [M. truncatula]PlastSol−0.160.040.24
8EGTA-1490gi|357512271Ribose-5-phosphate isomer. [M. truncatula]PlastSol-0.180.020.19
8EGTA-1794gi|157973737PSI reaction center IV A [Cicer arietinum]PlastSol−0.010.100.30
8LiCl-1156gi|537317peroxidase [M. sativa]ExtracSol−0.070.050.24
8LiCl-1163gi|537317peroxidase [M. sativa]ExtracSol0.030.130.22
8LiCl-300gi|657374433polygalacturonase non-catalytic ExtracSol0.000.060.20
8LiCl-340gi|537317peroxidase [M. sativa]ExtracSol−0.06−0.010.20
9CaCl2-1292gi|388508134Nod factor-bind. lectin [M. truncatula]GolgiMemb−0.22−0.150.14
9CaCl2-1940gi|388502358Xylogluc. EndotransglucosylaseExtracSol−0.31−0.010.11
9CaCl2-1945gi|357508519Xylogluc. EndotransglucosylaseExtracSol−0.100.010.21
9CaCl2-2042gi|5777611cysteine protease [M. sativa]ExtracSol−0.02−0.080.18
9CaCl2-2511gi|388502800Germin-like protein [Pisum sativum]ExtracSol−0.22−0.020.03
9EGTA-1376gi|357513539Stem 28 kDa glycoprotein [M. truncatula]ExtracSol−0.09−0.070.18
9EGTA-1391gi|357513539Stem 28 kDa glycoprotein [M. truncatula]ExtracSol−0.17−0.020.30
9EGTA-1661gi|357476945Acid beta-fructofuranosidaseLysoSol−0.150.010.18
9EGTA-956gi|537313plant/F18G18-200 protein [M. truncatula]ExtracSol−0.31−0.090.24
9EGTA-965gi|537313plant/F18G18-200 protein [M. truncatula]ExtracSol−0.23−0.050.24
10CaCl2-1586gi|537313plant/F18G18-200 protein [M. truncatula]ExtracSol−0.38−0.140.37
10CaCl2-1609gi|537313plant/F18G18-200 protein [M. truncatula]ExtracSol−0.45−0.200.42
10CaCl2-2570gi|388509752unknown [Lotus japonicus] (pfam03018)ExtracSol−0.16−0.130.26

Share and Cite

MDPI and ACS Style

Sergeant, K.; Printz, B.; Guerriero, G.; Renaut, J.; Lutts, S.; Hausman, J.-F. The Dynamics of the Cell Wall Proteome of Developing Alfalfa Stems. Biology 2019, 8, 60. https://doi.org/10.3390/biology8030060

AMA Style

Sergeant K, Printz B, Guerriero G, Renaut J, Lutts S, Hausman J-F. The Dynamics of the Cell Wall Proteome of Developing Alfalfa Stems. Biology. 2019; 8(3):60. https://doi.org/10.3390/biology8030060

Chicago/Turabian Style

Sergeant, Kjell, Bruno Printz, Gea Guerriero, Jenny Renaut, Stanley Lutts, and Jean-Francois Hausman. 2019. "The Dynamics of the Cell Wall Proteome of Developing Alfalfa Stems" Biology 8, no. 3: 60. https://doi.org/10.3390/biology8030060

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop